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I.  INTRODUCTION 


There  is  an  increasing  number  of  design  problems  vhich  involve  the 
coupling  of  magnetic  fields,  induced  currents  and  mechanical  behavior.  More¬ 
over,  many  of  these  problems  also  involve  thermal  effects  vhich  may  or  may 
not  be  coupled  with  the  magnetomechanical  behavior.  Among  the  types  of  machinery 
and  devices  for  vhich  such  coupled  behavior  is  significant  are  superconducting 
motors  and  generators  [l]  fusion  reactors  [2],  magnetic  levitation  transport 
systems  [2],  transmission  lines,  magnetomotive  tools  [3],  [U],  inductive 
storage  units,  MHD  devices  [2],  and  various  other  electromagnetic  devices  and 
svitches  [5].  Other  objects  subject  to  these  phenomena  are  structures  or 
devices  called  on  to  vithstand  lightning  or  electrical  shorting  in  high  voltage 
svitches  [5].  A  reviev  of  many  of  these  coupled  problems  in  magneto-solid 
mechanics  has  been  presented  by  Moon  [ 11 ] . 

The  computational  problems  involved  in  analysis  of  the  coupled  behavior 
of  such  structures  are  formidable.  Therefore,  prevailing  practice  has  been  to 
undertake  separate  analyses  of  the  magnetic,  thermal  and  mechanical  effects. 

A  magnetic  field  and  induced  current  analysis  is  first  undertaken  to  predict 
the  magnetic  forces  on  the  assumed  rigid  structure.  Temperature  fields  are 
next  predicted.  As  a  final  step,  analysis  for  the  uncoupled  mechanical  response 
is  performed. 

This  uncoupled  approach  may  be  inadequate  in  severed,  regards.  First, 
the  deformation  of  the  structure  will  alter  the  induced  currents  and  the 
concommitant  forces,  and  this  effect  vill  lead  to  altered  structural  response 
and  may  produce  structural  instabilities.  (For  example,  such  magnetomechanical 
instabilities  have  been  found  in  structures  for  proposed  fusion  reactors  [6], 
[7]).  Second,  in  dynamic  problems  the  eddy  currents  may  give  rise  to  thermo¬ 
elastic  effects.  Finally,  in  superconducting  problems,  the  vast  thermal 
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differences  and  gradients  will  have  a  pronounced  effect  not  only  on  the 
magnetic  behavior  but  also  on  the  material  nonlinearities  entering  the 
mechanical  behavior. 

The  goal  of  ongoing  research  at  Cornell  is  to  develop  numerical  methods 
for  the  coupled  analysis  of  forces,  currents,  and  stresses  in  thin  elastic 
structures  exposed  to  time -dependent  magnetic  fields.  In  particular  the 
project  is  directed  toward  time  dependent  magnetic  forces  and  motions  of 
elastic  plates  and  shells.  Both  applied  mathematical  analysis  and  finite 
element  numerical  computer  techniques  are  being  used.  Although  some  initial 
attempts  have  already  been  made  by  others  (Oden  and  Kelley  [8];  Miya  et  al. 

[9]»  Becker  and  Pillabury  [10]),  many  fundamental  difficulties  remain  to  be 
overcome  if  practical  analysis  tools  are  eventually  to  result.  Verification 
of  the  newly  developed  computer  codes  is  carried  out  by  direct  comparison  of 
the  calculated  results  with  experimental  data  and  is  carried  out  at  each 
stage  of  the  development.  The  experiments  employ  a  new  infrared  scanning 
technique  to  visualize  the  induced  eddy  current  patterns  in  the  structures. 

This  technical  report  is  an  account  of  the  first  year's  accomplishments  in 
this  effort. 

Coupled  magnetomechanical  problems  are  unique  in  that  the  electro¬ 
magnetic  field  equations  apply  both  inside  and  outside  the  structural  material. 
Thus  most  techniques  for  calculating  stresses  due  to  electromagnetic  fields 
require  a  three  dimensional  treatment  even  if  the  structure  is  one  or  two- 
dimensional  such  as  a  beam,  plate,  or  shell.  In  this  research  a  stream  function 
for  the  induced  current  in  the  conducting  structure  has  been  used  which  re¬ 
duces  the  magnetics  problem  to  a  two-dimensional  one  and  allows  one  to  use 
the  same  set  of  finite  elements  for  both  the  magnetic  and  elastic  deformation 
fields.  The  use  of  a  current  potential  or  stream  function  is  similar  to  vortex 
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mechanics  in  fluid  mechanics,  and  results  in  an  integro-differential  equation 
for  the  induced  currents  in  the  structure.  Once  these  currents  are  calculated, 
the  magnetic  forces,  stresses,  and  motions  of  the  structure  can  he  found. 

To  solve  the  integro-differential  equation,  a  finite  element  Galerkin  scheme 
has  been  developed  which  leads  to  an  algebraic  set  of  equations  for  the  induced 
currents.  The  resulting  global  matrix  is  complex,  nonsymmetric ,  and  fully 
populated.  The  usual  limited  connectivity  (banded  matrix)  finite  element 
representation  is  lost  in  contrast  to  conventional  structural  finite  element 
problems . 

To  date,  the  research  has  focused  on  the  response  of  flat  plates  to 
steady  state  harmonic  and  pulsed  magnetic  fields.  A  finite  element  code  for 
one-dimensional  magnetic  field  problems  has  been  successfully  developed 
which  can  calculate  the  induced  current  and  magnetic  field  in  the  plate,  the 
magnetic  pressure  distribution  on  the  plate,  and  the  induced  temperature  distri¬ 
bution  due  to  Joule  heating  in  the  plate.  The  computer  code  has  been  compared 
directly  with  experimental  results.  Comparison  of  calculated  induced  temperature 
fields  with  infrared  scanning  measurements  show  the  difference  to  be  within 
observational  error. 

A  two-dimensional  computer  code  for  the  induced  currents,  magnetic  field 
and  temperature  in  a  conducting  plate  has  also  been  written.  The  flow  of  eddy 
currents  around  a  crack  in  a  plate  is  one  problem  investigated  with  this 
program.  This  type  of  analysis  has  application  to  the  development  of  non¬ 
destructive  testing  techniques. 

In  the  next  chapter,  the  formulations  of  the  problem  are  summarized. 

This  summary  includes  the  basic  theory  of  coupled  magnetomechanics,  the  stream- 
function  method  for  eddy  currents  in  flat  plates  and  cylindrical  shells,  and 
the  coupled  magnetomechanic  formulation  for  flat  plates.  Chapters  III  and  IV 
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contain,  respectively,  descriptions  of  the  one-  and  two-dimensional  computer 
programs  for  eddy  currents  in  flat  plates.  Numerical  results  and  experimental 
verification  for  both  programs  are  presented  in  Chapter  V.  Finally,  conclusions 
based  on  developments  to  date  are  presented  in  Chapter  VI. 


I 


II.  FORMULATIONS 

In  this  chapter,  the  basic  theory  of  magnetothermomechanics  is  briefly 
reviewed.  The  formulation  of  the  eddy  currents  by  the  stream-function  method 
adopted  for  use  in  this  research  for  flat  plates  and  cylindrical  shells  is 
next  presented.  Finally,  the  coupled  magnetomechanical  formulation  for  flat 
plates  is  summarized. 

Review  of  Basic  Theory 

j 

Magnetothermomechanics  encompasses  a  rather  complicated  three-field 
problem.  The  coupling  of  the  three  fields  is  shown  schematically  in  Figure 
II-l  and  occurs  through  the  field  equations,  constitutive  equations,  and 
boundary  conditions.  A  number  of  formulations  of  the  governing  equations  and 
boundary  conditions  for  various  physical  specializations  have  been  presented 
[ll],  [12].  For  the  descriptive  purposes  of  this  report,  only  a  single  version 
corresponding  to  Figure  II-l  is  presented  here. 

When  a  deformable  electric  conductor  is  placed  in  an  electromagnetic 
field  B,  electric  currents  flow  in  the  solid  and  electric  and  magnetic  polari¬ 
zation  of  the  material  may  be  induced.  In  this  study  the  electric  polarization 
is  not  considered.  The  electric  currents  J  (sometimes  called  eddy  currents) 
and  magnetization  M,  can  then  interact  with  the  magnetic  field  to  produce  a 
body  force  distribution  on  the  solid  [13].  Thus  to  determine  the  stresses 
and  concommitant  deformation  of  the  structure,  the  distribution  of  electric 
currents,  magnetic  field,  and  magnetization  in  the  body  must  be  determined. 

In  this  research  only  nonferromagnetic  materials  (M  *  0)  such  as  aluminum  or 

copper  have  been  treated  to  date.  j 

In  addition  to  magnetic  forces,  the  flow  of  current  produces  heating 
in  the  solid  proportional  to  J  •  J  called  Joule  heating.  If  the  currents 
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Lorentz  body  force  (J  x  B) 


Figure  II-l.  A  Schematic  of  One  Form  of  Magnetothermomechanical  Coupling 
(E  =  electric  field ,  J  3  electric  current,  B  *  magnetic  flux, 

V  =  velocity  field,  M  =  magnetic  permeability,  a  -  coefficient 
of  thermal  expansion,  AT  =  change  in  temperature  from  refer¬ 
ence  temperature  ,  H  =  magnetic  intensity. )  , 
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are  confined  to  one  region  of  the  structure  then  differential  heating  can 
occur  and  thermoelastic  stresses  can  result. 

The  following  equations  and  boundary  conditions  apply  to  the  problem 
of  dynamic  thermoelasticity  coupled  to  the  quasistatic  electromagnetics  of  a 
medium.  The  medium  is  assumed  to  be  well  conducting,  non— polarizable  and 
non-magnet izable.  The  influence  of  strain  rate  on  temperature  distribution 
is  assumed  to  be  small  as  is  the  effect  of  temperature  on  material  properties. 
All  configurations  are  referred  to  a  Cartesian  system  (x^  x2>  x3)  fixed  in 
space.  The  assumption  of  quasistatic  electromagnetics  in  association  with 
dynamic  mechanical  behavior  is  Justified  by  the  fact  that  the  frequencies  of 
mechanical  vibrations  and  waves  are  generally  much  smaller  than  those  of 
electromagnetic  waves  of  the  same  wave  length. 

The  quasistatic  magnetic  field  equations  for  this  case  are  obtained 
from  the  Maxwell  equations  and  are  [13], 


7  x  E  »  -3B/3t 

(Il-la) 

7  •  B  -  0 
••  **■ 

(Il-lb) 

V  x  H  *  J 

(Il-lc) 

in  which  E  is  the  electric  field  density,  B  is  the  magnetic  flux  density, 

H  is  the  magnetic  intensity,  and  J  is  the  electric  current  density.  The 
electromagnetic  constitutive  equations  for  a  linear,  isotropic,  homogeneous, 
resistive,  nonferromagnetic*  material  are 


B  - 

•m 

li  H 

(ll-2a) 

J  » 

0[E  +  V  x  B] 

(Il-2b) 

*For  linear  ferromagnetic  material  one  replaces  by  VQVr,  y  >  1. 


-10- 


in  which  yQ  is  the  magnetic  permeability  of  vacuum,  a  is  the  electric  conduc¬ 
tivity  and  V  is  the  velocity  field.  For  a  zero  resistive  or  superconducting 
material,  Ohm's  lav  (II-2b)  is  replaced  by  an  equation  of  the  form 

E  *  aJ 

■«*» 

The  electromagnetic  boundary  conditions  are 

nx  I  E  +  V  x  b1  *  0  (Il-3a) 

**  l  —  **■’  *“Ji 

n  x  JhJ  -  0  (II-3b) 

n  x  [b]  *  0  (Il-3c) 


in  which  n  is  a  unit  normal  vector  to  the  boundary,  and  JJ¥j|  denotes  the 
jump  of  the  quantity  ¥  across  the  interface. 

The  mechanical  field  equations  are  the  equilibrium  equations 


3%l  3Ti1 

p  rr  -  sr  +  Fi ♦  <£  *  5>i 


(II-U) 


in  which  u  is  the  displacement  field,  is  the  stress  tensor,  F  is  the 
mechanical  body  force  intensity,  p  is  the  mass  density.  Equation  II-4  includes 
magnetomechanic al  coupling*  through  the  Lorentz  body  force  term  (J  x  B).  The 
mechanical  constitute  equations  for  a  linear,  isotropic,  nonferromagnetic 
material  is  given  by  the  Duhamel-Neuaann  lav. 


Tij  *  20e1^  *  a  “  8  (T  -  To)]  Sij 


(II— 5 ) 


*For  ferromagnetic  bodies  an  additional  force  must  be  added. 
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in  which  G  and  X  are  the  Lam£  parameters,  e^  in  the  strain  tensor,  T  is  the 
temperature,  Tq  is  the  initial  uniform  reference  temperature ,  0  ■  (3d  ♦  2G)a, 
and  a  is  the  coefficient  of  thermal  expansion.  Finally,  the  mechanical  boundary- 


conditions  are 


Tu  nj  '  fi  *  ffu  Eiaj  -  \  HA  siJnj °n  3t 


(II-6a) 


(II-6b) 


in  which  S  ■  ♦  S^,  and  and  u^  are  specified  surface  tractions  and 

surface  displacements. 

For  the  third  field,  the  temperature,  the  field  equations  are 


or  H  «•  32T 

^vJt  "  K  3xi  3x^ 


J  •  E 


(II-7) 


in  which  is  the  specific  heat,  and  K  is  the  thermal  conductivity.  The 
Joule  heating  term,  J  •  E,  provides  magnetothermal  coupling.  The  strain  rate 

•m  mm 

coupling  term  has  been  dropped  in  (II-7).  In  the  absence  of  heat  input 
euros s  boundaries  and  if  radiation  can  be  neglected  the  thermal  boundary- 


conditions  arv 


T  on  S. 


0  on  S„ 


(Il-8a) 
( II-8b ) 


in  which  S  *  S1  +  Sg,  and  T  is  the  specified  surface  temperatures. 

Equations  II-l  through  II-8  constitute  an  example  of  a  coupled  problem  of  the 
form  shown  schematically  on  Fig.  II-l.  A  practice  frequently  followed  is  an 
uncoupled  approach  in  which  first  the  electromagnetic  problem  is  solved  for 
"rigid"  conductors,  second  the  temperature  distribution  due  to  Joule  heating 
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is  ascertained,  and  finally  the  mechanical  problem  for  the  conductors  is  con¬ 
sidered  driven  by  the  known  temperatures  and  Lorentz  forces.  This  simplified 
problem  is  obtained  in  Figure  IX-1  by  omitting  the  dashed  coupling  lines  from 
the  schematic. 

For  dynamic  magnetic  fields  even  the  uncoupled  electromagnetic  problem 
is  difficult  to  solve  since  the  induced  currents  in  the  conductor  are  unknown. 
When  an  external  magnetic  field  B°  induces  eddy  currents  in  a  conductor  these 
currents  generate  their  own  self  magnetic  field  (see  Figure  II-2)  which  in 
turn  affects  the  eddy  current  problem.  Several  problems  of  this  nature  have 
been  solved  using  finite  element  methods  ([lb]  -  [17]).  These  approaches 
however  use  a  finite  element  discretization  both  inside  and  outside  the  conductor 
vhich  makes  the  extension  of  these  numerical  schemes  to  coupled  mechanical 
problems  difficult. 

Two  different  solution  approaches  to  the  eddy  current  problem  have  been 
undertaken.  The  first,  or  local  approach,  seeks  to  solve  the  differential  field 
equations  and  satisfy  all  the  boundary  and  initial  conditions  directly.  The 
domain  of  the  solution  is  the  entire  space  in  vhich  the  fields  act  and  includes 
the  current-source  conductors,  the  induced-current  conductors,  and  the  sur¬ 
rounding  space  (theoretically  out  to  an  infinite  distance).  The  second,  or 
nonlocal  approach,  uses  Integral  restatements  of  the  formulation  such  that 
only  the  bodies  of  the  conductors  themselves  need  be  considered  in  the  problem. 
Thus  this  approach  avoids  treatment  of  free  space  extending  to  infinity  and  may 
be  easier  to  extend  to  coupled  magneto-mechanical  problems. 

Briefly,  the  local  approach  for  the  uncoupled  problem  of  Eqs.  II-l  through 
II-8  above  is  most  conveniently  expressed  in  terms  of  the  magnetic  vector 
potential  A,  defined  by 


B  ■  V  x  A  (since  7  •  B  ■  0) 


(II-9) 


The  governing  equations  for  the  potential  A  in  the  conductor  become 

*** 

*  1  # 

0  A  +  -  7xVxA-Oux7xA  -  Jq  (II-10) 

in  which  J  is  the  initial  current  density  carried  by  the  conductor. 

Outside  the  conductors  O  ■  u  ■  0  and  A  must  satisfy 

*•» 

V  x  V  x  A  ■  0 

for  non-electromagnet ic  wave  problems.  Thus  the  solution  for  A  must  be 
pieced  together  for  regions  outside  and  inside  the  conductors  using  appropriate 
boundary  conditions.  It  is  difficult  to  couple  a  structural  finite  element 
scheme  for  the  deformation  of  the  conductors  using  this  local  formulation. 

A  Galerkin  finite  element  formulation  can  be  used  for  this  form  of  the 
problem  ( [lU]  -  [17]),  but  details  will  not  be  repeated  here. 

Alternatively,  a  nonlocal  formulation  for  the  eddy  current  problem  may 
be  constructed  by  the  integral  equation  method.  Here  the  vector  potential 
is  separated  into  two  independent  parts ,  A°  and  A1 : 

A  *  A°  +  A1  (11-11) 

corresponding  to  the  driving  and  induced  fields,  respectively.  Ohm's  lav 
for  the  rigid,  induced-current  conductor  may  be  written  as 

J1  ■  OE1  »  -0  ~  0  It  k°  ~  °  -  (11-12) 

where  J°  has  been  assumed  to  be  zero,  i.e.,  no  external  voltage  is  applied 
to  the  induced-current  conductor.  Within  the  conductor,  the  external  field 
satisfies  7  x  H°  »  0;  therefore,  A0  is  harmonic,  i.e.,  7^A°  *  0,  while  the 
induced  field  satisfies  the  vector  Poisson  equation 


because  7  x  *  A  If  retardation  is  neglected  due  to  sufficiently 
slow  variation  of  all  electromagnetic  quantities  with  time,  the  solution 
of  Eq.  (11-13)  for  A1  is 

-1  "  iff  /  \x  -Sx'|  d<Vol'>  (n-uo 

▼ 

which  gives  an  integral  equation  in  terms  of  J* 

i1  +$  f  ]r  -~r'T  Jt  i1  (£f)  d(Vol,)  “  -a  it  ?  (II-!5) 

v 

or  in  terms  of  A1 

i1  ♦  ft  /  T TT^T  *  aI  (r,)  a(Vol,)  -  -  (II-16) 


A  finite-element-like  treatment  of  this  formulation  of  the  eddy  current 
problem  has  been  presented  by  Sylvester  et,  al .  [l8]. 

Eddy  Currents,  Pressures  and  Temperatures  by  the  Stream  Function  Method 

When  conductors  are  thin-sheet  structures  such  as  plates  or  shells 
which  are  used  in  many  practical  configurations,  a  different  nonlocal  approach 
for  the  eddy  current  problem  suggested  by  Moon  [11 ]  can  be  of  greater  utility 
than  the  integral  equation  method.  This  involves  the  introduction  of  a  stream 
function  for  the  current.  If  the  current  is  of  relatively  low  frequency, 
the  induced  current  distribution  can  be  assumed  to  be  uniform  through  the 
thickness,  and  the  current  density  is  denoted  by  I  *  hJ  where  h  is  the  thick¬ 
ness  of  the  sheet.  Let  (o, 3)  be  a  two-dimensional  orthogonal  coordinate  system 
on  the  mid-surface  of  the  sheet  conductor  structure,  and  let  (x,y,z)  be  a  Cartesian 
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coordinate  system,  so  that 

x  ■  x(a,8).  y  *  y(a,B),  z  ■  z(a,B) 

The  infinitesimal  line  element  on  the  surface  is  given  by 

ds  *  f  da  +  g  dp 
in  which 

r2  ■  <|f>2  ♦  <|£>2  *  <|f>2 
*2  -  #2  *  #2  -  #2 

The  orthogonality  of  (a, 8)  coordinate  system  requires  that 

/3xw3xa  (^L\(^SL\  +  «  0 

3a  3?  3a3B  l3an3BJ 

Let  e  and  ea  be  unit  tangent  vector  along  a  *  constant  and  B  9 
—a  ~p 


respectively. 


«a 

**» 


1  3x  \ 
f  3a  \»  26 


1  &L 

f  3a 

1  3z  / 
f  3a 


1  3x  . 
g  »6  \ 

1 

g  38 

1  3z  / 
g  3K 


The  current  density  on  the  surface  can  be  written  in  the  form 

I  -  l“(a,B)2a  +  I0(o’6)  !b 

The  continuity  of  I  requires  that 

div  J  *  ffj  +  ^flS)1  *  ° 

This  suggests  the  use  of  a  stream  function  V  for  the  current 
This  function,  when  defined  as 


.a  1  3i  t8  _  1  3i 

1  *  g  38  ’  "  f  3a 


(II-17) 

(11-18) 

(II-19) 

(11-20) 
constant , 

(11-21) 

(11-22) 

(11-23) 


(11-24) 


identically  satisfies  the  continuity  condition.  The  current  flovs  along  lines 
of  ¥  *  constant.  The  time-varying  magnetic  field  is  the  sum  of  the  external 
field,  B°,  and  the  induced  field,  B^,  Therefore,  from  Faraday's  lav. 

Equation  II-l(a),  curl  E  =  -  ^  B,  one  may  write  a  governing  equation  for 
the  determination  of  ¥: 

curl  •  curl  (4m)  *  ah  -rr  (B°  +  B1) 

3t  -  -  (11-25) 

Only  the  normal  (n)  component  of  the  vector  equation  will  be  used 
because  the  in-plane  current  is  of  primary  interest.  The  induced  field  may 
be  obtained  from  the  Biot-Savart  law.  The  specific  derivations  are  given 
below  for  flat  plates  and  thin  cylindrical  shells.  Once  the  relationship 
between  the  induced  field  and  the  stream  function  is  obtained,  ¥  may  be 
determined  from  Eq.  ( 11-25).  The  current  can  be  recovered  from  Eq.  (II-2U), 
and  magnetic  pressure  and  temperature  then  obtained  for  the  calculated 
current. 

Flat  plate.  The  derivation  for  the  case  of  a  flat  plate  has  been 
given  by  Moon  [ill.  The  details  are  not  repeated  here;  only  the  essential 
relationships  are  recapitulated. 

If  the  mid-plane  of  the  plate  is  chosen  as  the  xy  plane,  a*x, 

1  8*y,  and  one  obtains 

I  *  V  x  (tpn)  *  -n  x  VW 

j  ~  -  (11-26) 

I  and 

V2^  *  Oh  ^  (B®  +  B*)  (II-27) 

For  the  determination  of  E*,  the  current  stream  function  ¥  is  assumed 

z 

to  be  defined  throughout  three-dimensional  space,  and  constant  across  the 

.  ■— ~r 

I - -  — - .  .  . . . . . . 


sheet  thickness.  Application  of  the  Biot-Savart  law  and  the  divergence 

theorem  yields  B*  for  the  midplane  of  the  plate  (z  =  0) : 
z 


BI(x,y,o)  -  -  fa  \ 

area 


i|,(x',y')  a(area') - (II_28, 

t(x  -  i1)2  ♦  (y  -  y')2  ♦  (i?)2]3'2 


Y  has  been  set  to  zero  on  the  boundary  of  the  plate  in  the  above 

derivation.  When  the  plate  is  multiply  connected,  different  constant 

values  of  Y  will  be  assigned  on  each  interior  boundary.  These  boundary 

conditions,  together  with  the  known  B°  and  Eq.  (11-28),  are  then  used  in 

z 

( 11-27)  for  the  determination  of  the  stream  function  Y.  Once  Y  is  deter¬ 
mined,  the  current,  pressure,  and  temperature  may  be  calculated  from  ( 11-26) 
and  the  following  two  equations 


P  ■  I  x  (B°  +  B  ) 


(I-I)dt 


(11-29) 


(11-30) 


in  which  P  is  the  pressure,  T  the  temperature  rise,  p  the  density,  and 
the  specific  heat  of  the  conductor. 

Cylindrical  shell.  The  stream  function  method  has  been  extended  to 
the  geometry  of  a  cylindrical  shell.  In  this  case  the  coordinate  system 
is  chosen  as  shown  in  Fig.  II-3.  For  a  fixed  radius  p  *  a,  one  has  a  =  <}>, 
8  =  z.  The  stream  function  Y  may  be  introduced  such  that 


I-Vx  (i|»p)  $, 


( ii-3i : 


t  a  3i  i  =  _  i  M 

<p  3z  ’  z  P  ^ 

The  normal  component  of  ( 11-25)  in  this  case  reads 

±  it .  i± . «,  ♦  b2) 


( n-32 ; 


(11-33: 
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Figure  II-3.  Notation,  Coordinates, 
Cylindrical  Shell. 
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To  calculate  Bx  one  starts  from  the  Biot-Savart  law 
P 


?I(r>  -  *  J 


I(r')  *  (r  -  r') 


d(Vol. ' ) 


volume 


e  =  <j>  -  <f>* 


(11-3U) 


(n-35) 


R  *  (pcos0  -  p')  e  ,  +  psin9  efl,  +  (z  -  z')  e  , 

**  P  *0  *vZ 


Using  Eq.  ( 11-32)  and  the  relationship 


(11-36) 


7iW3]rT 


(11-37) 


one  has' 


bp(s)  =  [(Hr)(ib'iRf)  +  itf)d(vo1-' 

vol. 

The  integrand  can  he  written  in  the  form 

V  V1  lip  {W  1bT) 


(11-38) 


(11-39) 


The  second  term  vanishes  because  =  0.  Using  an  identity  to  replace 
the  first  term,  one  obtains  the  expression  for  B^  as 


Bfc>  “  sik  f  <Y (*vi|ir>  -  *7i  ft1  d<TOl- 1 

vol. 

7X  • iir> 1(TOl-> 


( 11-1*0) 


If  one  vises  the  divergence  theorem  on  the  second  integral  and  sets  ¥  =  0 
on  the  end  edges  of  the  cylinder,  one  has 


■  h  *  ♦  rar  \  *  ^  ur  i(are‘'  ’ 


outer 


-H-  4;  -r  -i—  d(area* ) 

4TTh  y  3p*  |Rj 

inner 


(II-Ul) 


I 


On  the  middle  surface  of  the  cylinder  one  obtains,  by  setting  p =  a 


in  the  above  expression. 


Bp(a’*’z)  “h^sSh  j 


^(a,^'  ,z*  )[acos0  -  (a  +  jjO](a  +  jjOdOdz' 

.  {[acos9  -  (a  +  |)]2  +  a2sin20  +  (z  -  z')2}3^2 
outer  2 


-  nfe  \ 


u  1. 

tj)(a,<(i'  ,z'  )[acos9  -  (a  -  r-)](a  -  v-)d0dz’ 


inner 


{[acos9 


~t  hn2  2.2,  ^372  (II-1*2) 

-  (a  -  «)]  +  a  sin  9  +  (z  -  z' )  J 


when  I  is  independent  of  the  variable  z ,  the  self  field  becomes 

2ir  . 

t  ^(9)[acos9  -  (a  +  •%)] 

^(..*.0)  -  s* -  <t- *  5-)  |  77777 w -  2a(a  7 fw de 

O  2  ^ 

2ir  h 

,  .  r  ij»(9)[acos9  -  (a  -  -r)} 

- § - - - d0  (II-U3) 

»  h  2  J  a2  ♦  (a  -  |)2  -  2a(a  -  |)coaS 


One  can  show  that  when  ¥  =  c,  a  constant  everywhere,  or  I  =  0, 


then 


2lT  [acoa0  -  (a  +  jjOKa  +  j) 


_I  uc  r.  .  1  f  f  LaC°30  -  lat?J;a+2' 

'o’*  &  W  J  {[acosd  -  (a  +  ;r)]2  +  a2sin20  +  "(z  -  a’)2}^2 


d9dz 


-oo  o 


2ir 


-M  | 


[acos9  -  (a  -  |-)](a  -  |-) 


{[acos0  -  (a  -  jj-)]2  +  a2sin20  +  (z  -  z')2}^2 
o  & 


d0dz 


^  (1  -  1  +  0) 


=  0 

as  it  should  be. 

The  same  boundary  conditions  as  in  the  flat  plate  case  also  apply  to 
cylindrical  shells.  The  current  may  be  calculated  from  ( 11-32)  after  ¥  is 
determined.  Pressure  and  temperature  may  be  computed  using  (11-29)  and  (11-30). 


Coupled  Magnetomechanics 

The  complete  set  of  field  equations  and  boundary  conditions  for 
coupled  magnetothermomechanics  has  been  presented  in  the  first  section. 

The  finite  element  Galerkin  technique  may  be  applied  to  reduce  them  to  a 
set  of  nonlinear  algebraic  equations.  Although  all  the  coupling  effects 
are  preserved  in  the  finite  element  equations,  this  approach  does  not  seem 
to  be  practically  feasible.  The  reasons  are  obvious.  Firstly,  the  three- 
dimensional  nature  of  the  problem  cannot  be  avoided.  Both  the  conductor 
and  the  surrounding  space  need  to  be  treated,  and  some  fictitious  boundary 
conditions  assumed  at  some  fair  away  distance.  This  makes  the  storage  and 
computational  efforts  needed  to  solve  these  nonlinear  equations  formidable. 
Secondly,  the  geometric  and  time  characteristics  of  the  problem  cannot  be 
easily  represented.  Consequently,  it  will  be  difficult  to  study  the  roles 
of  some  physical  parameters  in  practical  problems.  The  complete  field 
approach  is  not  emphasized  in  this  research  work;  therefore,  the  finite 
element  formulation  of  this  approach  will  not  be  presented  here. 

The  stream  function  method,  on  the  other  hand,  averages  the  current 
through  the  thickness  and  represents  it  on  the  mid-surface  of  the  sheet-type 
conductors.  The  problem  is  therefore  reduced  to  a  two-dimensional  one.  This 
method  employs  the  current-sheet  nature  of  the  eddy-current  distribution, 
and  is  parallel  in  spirit  to  the  plate  and  shell  theories  of  elasticity. 

A  definition  of  the  mid-surface  is  needed  in  this  method.  All  the  physical 
quantities  are  calculated  and  represented  on  it.  This  nature  of  the  stream 
function  approach  makes  it  intrinsically  Lagrangian.  The  incorporation  of 
the  large  deformation  analysis  of  structure  and  the  transient  eddy  current 
calculation  seems  natural  in  this  light,  although  some  kind  of  surface  inter¬ 


polation  may  be  needed  to  represent  the  deformed  structure  in  the  intermediate 


stages.  A  possible  solution  procedure  for  large  deformation  analysis  of 
the  coupled  magnetomechanics  using  stream  function  method  is  shown  in 
Fig.  II-U. 

From  the  practical  point  of  view,  magnetomechanical  problems  may  be 
classified  into  several  groups  according  to  the  difference  in  time  scales 
between  mechanical  and  electromagnetic  field  problems.  In  the  impulse 
loading  problems,  such  as  electromagnetic  metal  forming,  impulse  durations 
of  a.  io1  ysec  and  strain  rates  of  ^  lO^sec  ^  have  been  reported  [19  ] .  Since 
the  fundamental  period  of  the  conductor  structure  is  much  higher  than  the 
duration  of  the  pulse,  these  problems  are  best  treated  as  initial  velocity 
problems.  The  nonlocal  eddy  current  effect  is  expected  to  be  important  in 
this  case  and  is  definitely  needed  for  the  accurate  calculation  of  the 
pressure.  The  mechanical  part  of  the  problem  in  this  case  is  usually  one 
of  dynamic  plasticity. 

When  the  duration  of  the  pulse,  or  the  period  of  the  harmonic  driving 
current,  is  comparable  with  the  fundamental  period  of  the  structure,  most  of 
the  coupling  will  occur  through  the  displacement  of  the  structure.  The  non¬ 
local  effect  will  be  small  in  this  case,  because  the  eddy  current  will  have 
enough  time  to  diffuse.  A  much  simpler  problem  may  then  be  obtained  by 
dropping  the  nonlocal  effect  of  the  eddy  current.  If  the  displacement  is 
small  so  that  the  original  configuration  may  be  used  for  reference  for  both 
structural  and  eddy  current  problems,  two  sets  of  linear  matrix  equations  may 
be  written  for  structural  and  eddy  current  problems  separately.  The  coupling 
of  the  magnetic  field  and  the  motion  of  the  structure  in  this  case  appears  on 
the  right  hand  side  of  the  equations  only. 

The  most  difficult  coupling  problem  comes  from  the  cases  in  between  the 
above  two  extremes,  i.e.,  when  the  frequency  of  the  field  is  high  and  the 
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duration  of  the  pulse  long  enough  to  cause  large  deformations.  Examples  of 
these  problems  may  he  found  in  the  design  of  fusion  reactors  during  the 
starting  and  shut-off  periods.  The  full  power  of  the  algorithm  shown  in 
Fig.  II-U  may  be  needed  for  these  problems. 

Other  problems  of  interest  include  the  coupling  between  field  and 
the  rigid  body  motion  of  the  conductor,  such  as  that  of  magnetic  transportation . 
A  Judicious  choice  of  the  type  of  magnetomechanical  coupling  is  therefore 
essential.  For  plate  and  shell  type  problems,  the  stream  function  method  is 
believed  to  provide  a  more  convenient  way  to  model  the  various  coupling 
effects. 
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III.  ONE-DIMENSIONAL  CODE  FOR  FLAT  PLATES 


In  this  chapter  the  finite  element  Galerkin  formulation  for  the 
problems  of  one-dimensional,  steady-state  harmonic  currents  in  long  plates 
is  described.  The  treatment  of  the  nonlocal  integral  term  is  briefly 
discussed.  The  capabilities  of  program  ONED,  a  fortran  code  developed 
based  on  this  formulation,  is  then  introduced.  The  input  and  output  of 
program  ONED  are  finally  described. 

Formulation 

When  the  current  I  is  independent  of  y,  eq.  ( 11-27)  has  the  following 


form: 


A 

♦  i-ouh  f - - — -  ax'  •  Oh 

v?  *  2*  J  (x'  -  x)2  *  (§>2 


-A 


3B° 

--r  (in-i) 

3t 


For  3teady-state  harmonic  currents  in  the  plate,  this  equation  may  be  non- 
dimensionalized  into 

2  ® 

-  1R$  +  licr  [  ■  -  %  ■ — t  *  1RB(x) 

d7  (c-x)2  +  } 


(III-2) 


in  which  x  and  5  are  distances  across  the  width  nondimens ionali zed  with  respect 


_  i  /lit.  _  ,  •  »  /  »  »  I  fii t . 

to  the  thickness  h;  ^  *  $I°e  ;  B°  *  (y/A)B(x)l°e  ;  the  magnetic  Reynolds 

z 


2  2 
number  R  *  w^ah  ,  which  is  related  to  the  skin  depth  6  through  R  *  2(h/6)  ; 


1°  is  the  current  in  a  nearby  exciter  coil;  2A  is  the  width  of  the  plate; 
and  2a  the  nondimensionalized  width. 

The  finite  element  Galerkin  method  is  used  to  solve  the  integro-differ- 
ential  eq.  (III-2).  $  is  approximated  globally  and  locally  by  piece-wise 


linear  models 
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in  which  G  is  the  total  number  of  nodal  points,  the  superscript  E  denotes 
the  E-th  element,  and  are  the  global  interpolation  functions  generated  from 
the  local  linear  element  shape  functions  N^.  The  linear  algebraic  equations 
for  each  element  are 
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in  which 


ik  ■ R  [ N!  i dx 

dC  ;  Rj  -  R  [  NjBdx 
E 


(III-6) 


the  weighting  functions  used  in  calculating  the  integral  term  is 


w'u) 
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_£(x) 


dx 


(HI-7) 


Both  (III-6)  and  (III-7)  are  integrated  analytically.  The  resulting  global 
matrix  is  complex,  nonsymmetric ,  and  fully  populated.  The  usual  limited 
connectivity  (banded  matrix)  finite  element  representation  is  lost  due  to 
the  nonlocal  nature  of  the  problem. 

When  the  frequency  of  the  external  magnetic  field  is  low,  the  effect 
of  the  self  field  is  small.  The  solution  in  this  case  will  approach  the 
local  solution  obtained  by  dropping  the  second  and  third  terms  on  the  left 
hand  side  of  (III-2).  When  the  frequencies  is  high,  the  distribution  of  the 
eddy  current  is  primarily  influenced  by  the  self  field.  The  solution  will 
approach  the  one  obtained  by  dropping  the  first  term  on  the  left  of  (III-2). 
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The  solution  so  obtained  is  an  asymptotic  value  for  the  high  frequency  limit, 
and  is  called  the  image  solution  in  the  literature. 

Capabilities  of  the  Program  ONED 

A  Fortran  program  ONED  has  been  developed  based  on  eqs.  (III-5)  -  (III-T). 
It  is  designed  to  perform  the  folloving  calculations: 

(1)  local  solutions  of  stream  function,  current,  temperature,  and 
pressure , 

(2)  nonlocal  solutions  of  stream  function,  current,  temperature, 
and  pressure, 

(3)  image  solution  of  stream  function,  current,  tenq>erature ,  and 
pressure. 

To  perform  the  calculations  necessary  to  construct  a  spectrum,  0NED  performs 
the  above  calculations  for  any  number  of  magnetic  Reynolds  numbers  in  one 
run.  Multiple  load  cases  (up  to  10)  are  allowed.  Each  load  case  will  be 
solved  for  all  the  magnetic  Reynolds  numbers  specified  for  that  run. 

Presently,  ONED  can  handle  uniform  magnetic  field  cases  and  cases 
where  the  fields  are  generated  by  any  number  of  conductor  wires  parallel  to 
the  plate.  The  strengths  and  senses  of  the  currents  and  the  positions  of  the 
wires  may  be  different.  Space  has  also  been  left  to.  include  other  types 
of  exciting  fields  in  the  program. 

Input 

Program  ONED  takes  the  geometric  information  of  the  plate  and  magnetic 
field  source  information  a3  input.  The  total  numbers  of  nodal  points,  load 
cases,  parameter  sets,  and  the  plotting  option  of  the  output  must  be  specified 
on  the  master  input  card.  Nodal  point  information  may  be  generated  for  portions 


of  the  mesh  that  are  uniform. 
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The  type  of  source  of  magnetic  field  is  specified.  No  other  information 
is  needed  for  uniform  magnetic  field.  For  the  parallel  wire  cases,  the 
positions,  senses,  and  the  relative  strengths  of  the  currents  in  each  wire 
must  he  input  subsequently. 

All  the  input  data  must  be  nondimensionalized  according  to  the  con¬ 
ventions  presented  in  the  formulation. 

Output , 

Values  of  the  stream  function,  eddy  current,  temperature  induced  in  a 
half-cycle  of  the  current,  and  time-averaged  magnetic  pressure  exerted  on  the 
plate  are  produced  as  output  of  the  program.  The  stream  function  and  eddy 
current  are  calculated  in  complex  form.  The  stream  function  values  are 
calculated  at  the  nodal  points  of  the  mesh,  while  the  current,  temperature , 
and  pressure  are  evaluated  at  the  centroids  of  the  elements.  The  stream 
function  and  eddy  current  are  given  in  complex  form,  and  the  modulus  and 
phase  angle  are  calculated  for  the  eddy  current.  Values  of  the  modulus  and 
phase  angle  for  different  frequencies  may  be  used  to  generate  the  spectrum 
of  the  current  at  any  point  on  the  plate.  These  curves  and  the  spectrum  of 
the  exciting  current  may  be  used  to  calculate  the  transient  current  at  the 
point  by  fast  Fourier  transform  techniques. 

All  the  outputs  of  the  program  may  be  plotted  out  using  a  printer-plotter 
subroutine.  This  capability  is  optional  and  can  be  specified  in  the  input. 

In  this  program  the  local  solutions  are  printed  out  and  plotted  parallel 
to  the  nonlocal  solutions  for  comparison.  When  there  are  several  load  cases, 
the  results  are  printed  out  in  groups  in  the  sequence  of  the  load  numbers. 

In  the  multiple  parameter  (frequency)  cases,  this  is  done  for  each  Reynolds 
number.  Title  lines  will  be  printed  for  each  parameter  to  distinguish  the 
different  groups  of  output.  Image  solutions,  which  are  independent  of  the 


IV.  TWO-DIMENSIONAL  CODE  FOR  FLAT  PLATES 


The  finite  element  Galerkin  formulation  for  the  two-dimensional  problems 
of  steady-state,  harmonic  currents  in  flat  plates  is  presented  first.  The 
treatment  of  the  nonlocal  integral  term  is  then  described  in  some  detail. 
Program  TWOD,  developed  based  on  this  formulation,  is  then  introduced.  The 
input  and  output  of  TWOD  are  described  at  the  end. 

Formulation 

For  steady  state,  harmonic  current  in  a  flat  plate  plate,  ( 11-27)  may 
be  nondimens ionali zed  into  the  following  form 


V^4>  -  i2irR4»  +  iR 


area 


[( x-  o  +  (y  -  n)  +  1] 


3/2 


d£dr|  =  i2irR9(x,y)  (IV-1) 


in  which  the  coordinates  are  nondimensionalized  with  respect  to  half  the 
thickness  ( j) ;  ip  =  ;  B°  *  (^)0e*Wt;  the  magnetic  Reynolds  number 

is  R  =  ,  which  is  related  to  skin  depth  6  through  R  ®  (y)2;  and  B 

is  the  reference  magnetic  field. 

The  finite  element  Galerkin  method  is  used  to  solve  (IV-1).  <f>  is 

approximated  globally  and  locally  by 
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in  which  G  is  the  total  number  of  nodal  points,  E  denotes  the  Eth  element, 

are  the  quadratic  global  interpolation  functions  generated  from  the  local 

E 

element  shape  functions  N^.  Six-node  triangular  elements  are  used  here. 

The  local  element  shape  functions  are  all  quadratic  in  this  case.  The 
element  algebraic  equations  are 
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in  which 


4  ■  A  f  ^  “* 

ae 

*5,  - «  K  <  “* 

ae 

Qjk  3  -R  |{  V*»n)  (S,n)d5dn 

area 

■  -2irR  |  N^9dAE 


The  weighting  function  is  given  "by 
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(IV-5) 


Because  of  the  numerical  difficulty  associated  with  the  sharp  variation 

of  the  Kemal  function  in  (IV-5),  the  weighting  function  is  calculated 

analytically  within  the  element  and  numerically  outside  the  element. 

Eq.  ( IV-Uc)  is  then  integrated  entirely  by  numerical  quadrature. 

E 

Since  Nj  is  quadratic,  six  basic  integrations  with  numerators 
2  2 

l,x,y,x  ,xy,  and  y  in  (IV-5)  need  to  be  performed.  To  simplify  the 
integration,  each  element  is  first  mapped  onto  a  standard  triangle,  which 
is  independent  of  the  coordinates  of  the  nodes  of  the  element.  The  inte¬ 
grations  are  then  performed.  The  expressions  obtained  are  then  used  for 
the  calculations  of  (IV-5)  for  all  the  elements. 

The  resulting  global  matrix  is  again  complex,  nonsymmetric ,  and  fully 
populated.  Local,  nonlocal,  and  image  solutions  can  be  calculated  just  as 
in  the  one-dimensional  case.  Eddy  currents  can  be  calculated  through  numerical 
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differentiation,  and  magnetic  pressure  and  temperature  evaluated  after 
the  stream  function  is  solved. 

Capabilities  of  the  Program  TWOD 

A  fortran  program  TWOD  has  been  developed  based  on  the  formulation 
of  eqs.  (IV-3)  -  (IV-5).  It  calculates  the  local  and  nonlocal  solutions  of 
stream  function,  eddy  current,  temperature,  and  pressure.  As  of  this  writing, 
the  image  solution  and  two-dimensional  graphic  output  capabilities  remain  to 
be  implemented.  Uniform  magnetic  field  and  fields  due  to  any  number  of 
magnetic  dipoles  can  be  handled.  Magnetic  fields  generated  from  some  types 
of  coils  of  interest  will  be  added.  Capability  of  analyzing  for  multiple 
frequencies  is  also  to  be  implemented. 

Input 

The  geometry  of  the  plate  and  the  description  of  the  external 
magnetic  field  are  the  two  basic  forms  of  data  needed  by  program  TW0D. 

The  total  numbers  of  nodal  points,  load  cases,  and  elements  need  to  be 
specified.  Coordinate  and  boundary  condition  must  be  given  for  each  input 
node.  Intermediate  nodes  may  be  generated  for  any  groups  of  nodal  points 
that  are  uniformly  spaced.  Element  information  may  also  be  generated.  Although 
only  six-node  triangles  are  included  in  the  present  version  of  TW0D,  the 
program  has  been  structured  so  that  other  types  of  elements  may  be  added. 

Element  group  information  and  the  master  card  for  each  element  group  there¬ 
fore  need  to  be  input  too. 

The  program  allows  for  different  orders  of  numerical  integrations. 

Six-  and  seven-point  formula  are  now  provided.  The  order  may  be  specified 
on  the  master  element  group  card. 

Presently  two  types  of  magnetic  fields  may  be  analyzed:  the  uniform 


field  and  the  magnetic  field  due  to  any  number  of  dipoles.  For  uniform 
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field  no  other  information  is  needed.  For  dipole  field  the  positions  and 
the  relative  strengths  and  senses  of  each  magnetic  dipole  must  he  given 
subsequently. 

Output 

Values  of  the  stream  function,  eddy  current,  temperature  induced  in  a 
half  cycle  of  the  current,  and  time-averaged  magnetic  pressure  are  produced 
as  output.  The  stream  function  is  calculated  at  the  nodes  of  the  finite  element 
mesh.  Current,  temperature,  and  pressure  are  evaluated  at  the  centroid  of  each 
element.  The  stream  function  and  current  are  calculated  in  complex  form.  The 
modulus  and  phase  angle  of  the  current  are  evaluated  in  the  interest  of  spectral 
analysis  for  the  calculation  of  transient  currents. 
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V.  NUMERICAL  RESULTS  AND  EXPERIMENTAL  VERIFICATION 

The  computer  programs  described  in  the  preceding  two  chapters  have 
been  applied  to  a  number  of  problems  to  test  and  verify  their  validity  and  to 
demonstrate  their  utility  for  problems  of  scientific  and  practical  interest. 

A  portion  of  the  verification  process  is  the  comparison  of  the  computed 
results  with  experimental  measurements.  This  chapter  consists  of  a  summary 
report  of  these  numerical  results. 

One-Dimensional  Analysis 

The  first  test  of  the  numerical  techniques  discussed  in  the  previous 
sections  was  the  calculation  of  induced  currents  in  a  long  rectangular  plate 
by  a  two-dimensional  magnetic  field  (Figure  V-l).  Near  the  center  of  the 
plate  the  currents  are  one  dimensional  or  are  parallel  to  the  long  edges  of 
the  plate.  The  external  magnetic  field  used  in  the  test  cases  was  that  due 
to  parallel  current  filaments  above  the  plate  and  a  uniform  time  dependent 
magnetic  field.  In  addition  to  the  induced  current  distribution  across  the 
plate,  the  induced  temperatures  and  magnetic  pressure  distributions  were 
calculated. 

The  following  objectives  were  met  with  the  one-dimensional  program: 

a)  Comparison  of  low  magnetic  Reynolds  number  (R)  results  with 
direct  quadrature  of  local  theory. 

b)  Comparison  of  high  R  results  with  the  image  method. 

c)  Demonstration  of  the  importance  of  nonlocal  theory  for  moderate 
frequencies  or  R. 

d)  Comparison  of  finite  element  calculations  with  experimental 
infrared  measurements. 

e)  Combination  of  the  finite  element  (FE)  results  for  different 


field  frequencies  with  the  fast  Fourier  transform  (FFT)  to 
predict  dynamic  currents  and  pressures. 
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Figure  V-l.  Induced  Currents  in  a  Long  Rectangular  Plate. 
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f)  Demonstration  of  the  importance  of  the  edge  effect  in 
increasing  the  current  density  and  magnetic  pressure. 

g)  Calculation  of  the  magnetic  forces  on  a  long  plate  due  to 
a  tilted  coil. 

h)  Examined  the  effect  of  a  handed  matrix  as  an  approximation 
to  the  full  matrix  or  complete  nonlocal  solution. 

Each  of  these  will  he  described  in  the  following  subsections. 

Comparison  of  FEM  with  Direct  Quadrature  for  R«l.  When  R«1 
one  may  drop  the  nonlocal  terms  in  eq.  (III-2)  and  obtain  the  local  theory 
for  the  stream  function,  i.e. 

V2<|)  *  iRB(x)  (V-l) 

The  solution  for  the  one-dimensional  case  can  be  found  by  direct  integration 
of  a  given  B(x).  Comparison  of  the  low  Reynolds  number  solution  for  the 
finite  element  method  (FEM)  and  direct  quadrature  is  shown  in  Figure  V-2  for 
a  pair  of  current  filaments  centered  above  a  long  conducting  plate.  The 
results  show  excellent  agreement. 

Comparison  with  the  Image  Method.  When  the  frequency  is  high,  i.e., 
R«l, the  nonlocal  effects  act  to  prevent  the  magnetic  field  from  penetrating 
the  plate.  In  this  limit  the  solution  can  be  approximated  by  considering  an 
image  coil  below  the  plate  with  serves  to  cancel  out  the  normal  component 
of  the  total  magnetic  field  {Figure  V-3).  The  results  of  the  one-dimensional 
FEM  with  the  image  coil  calculation  are  shown  in  Figure  V-l*  and  again  show 
excellent  agreement. 

Nonlocal  vs.  Local  Theory  for  Eddy  Currents.  A  comparison  of  the 
local  and  nonlocal  theories  for  eddy  current  distribution  is  shown  in 
Figure  V-5  for  different  R  or  frequencies.  It  is  clear  that  even  though 
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the  skin  depth  may  be  3 ever  ad  times  the  thickness  (R  *  0.1,  6  *  k.5h),  the 
nonlocal  distribution  differs  significantly  from  the  local  theory.  The 
importance  of  the  nonlocal  effects  is  not  often  recognized  as  evidenced 
by  the  number  of  eddy  current  calculations  which  neglect  the  nonlocal 
contribution  [  16]  [  20]  • 

A  more  fundamental  difference  between  the  nonlocal  and  local  effects 
concern  the  magnetic  pressure.  When  the  nonlocal  or  self  field  of  the 
induced  currents  is  neglected, the  average  force  on  the  plate  for  a  sinus¬ 
oidal  current  is  zero.  The  nonlocal  or  self- field  effects  are  necessary  to 
obtain  a  time-averaged  force  on  the  plate  between  the  sources  of  the  external 
field  (such  as  colls  or  current  filaments)  and  the  currents  in  the  plate. 
Comparison  of  FE  Calculations  with  Experimental  Measurements.  An 


important  feature  of  this  research  program  has  been  the  experimental  veri¬ 
fication  of  the  numerical  calculations.  To  check  the  calculations,  measure¬ 
ments  of  one  of  several  physical  quantities  along  the  surface  of  the  plate  must  be 
made  such  as  electric  or  magnetic  fields,  induced  temperature,  or  magnetic 
pressure.  Measurement  of  electric  or  magnetic  fields  associated  with  the 
induced  currents  in  the  plate  involve  taking  data  of  a  sufficient  number  of 
points  to  map  out  the  current  distribution.  However,  in  recent  years  an 
infrared  scanning  technique  has  been  developed  which  can  sense  the  incremental 
temperature  distribution.  This  technique  is  based  on  the  fact  that  for 
small  times  after  the  induced  currents  are  generated,  the  heat  conduction 
can  be  neglected  and  the  induced  temperature  is  proportional  to  the  square 


of  the  current  density.  The  heat  equation  is  given  by 

9t  O 


(V-2) 
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in  which  <  is  the  thermal  conductivity,  and  c  is  the  heat  capacity  of  the 
plate.  Under  appropriate  conditions  on  the  gradient  of  temperature 
and  time, one  may  neglect  the  first  term  and  write 

1  4  £  f  ^  (V-3) 

o 

To  simulate  a  one-dimensional  problem,  a  rectangular  multiturn 

induction  coil  was  placed  parallel  to  a  flat  aluminum  plate  and  induced 

temperatures  were  measured  across  the  plate  under  the  middle  of  the  coil 

(Figure  V-6).  The  two-dimensional  infrared  scanner  used  is  described  in 

2 

Appendix  A.  An  infrared  thermogram  of  the  temperature  or  J  pattern  in  the 
plate  due  to  a  pulsed  rectangular  coil  is  shown  in  Figure  V-7.  The  induction 
coil  has  damped  sinusoidal  currents.  Comparison  of  the  qualitative  behavior 
of  the  measured  temperature  and  the  calculated  temperature  distribution 
along  the  plate  using  an  effective  frequency  is  shown  in  Figures  V-8  and 
V-9  for  both  high  and  low  frequencies.  Both  calculated  and  measured  data 
have  been  normalized.  Excellent  qualitative  agreement  can  be  seen  for  the 
coil  centered  above  the  plate  (Figure  V-8)  and  for  the  coil  near  the  edge 
of  the  plate  (Figure  V-9). 

Quantitative  agreement  is  difficult  because  the  calculations  are  for 
a  harmonic  excitation  field  and  the  experiment  used  pulsed  current  in  the 
coil.  This  problem  was  solved  by  calculating  the  induced  current  density, 
as  a  function  of  frequency,  and  using  a  fast  Fourier  transform  (FFT)  to 
calculate  the  induced  current  as  a  function  of  time,  as  shown  in  Figure 
V-10.  This  data  was  then  integrated  using  equation  (V-3)  above  to  calculate 
the  temperature  as  a  function  of  time  at  the  point  of  maximum  temperature  in 
the  plate.  Using  the  maximum  temperature  in  time  along  with  the  calculated 
distribution  as  in  Figure  V-10  a  quantitative  comparison  of  calculated  and 


Figure  V-6.  Method  of  Detecting  Eddy  Current  Flow  Patterns  in 
Thin  Structures. 


Figure  V-7.  Photograph  of  an  Infrared  Thermogram  for  Currents  in  a 
Plate  as  Induced  "by  a  Pulsed  Rectangular  Coil. 
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Figure  V-8.  Comparison  of  FE  Solution  and  Measured  Temperatures  for  a 
Pair  of  Current  Filaments  Centered  Above  a  Long  Conducting 
Plate  (R  =  0.071). 
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Figure  V-9.  Comparison  of  FE  Solution  and  Measured  Temperatures 
for  a  Pair  of  Current  Filaments  Nearer  One  Edge  of  a 
Long  Conducting  Plate  (R  =  0.071*). 


measured  temperatures  were  made  as  shown  in  Figure  V-8.  As  might  be 
expected  the  measured  values  are  about  20/5  below  the  calculated  values.  This 
may  be  accounted  for  by  heat  transfer  effects,  and  the  fact  that  the  observed 
temperature  on  the  back  side  of  the  plate  may  be  different  than  the  mean 
temperature  across  the  plate. 

It  is  believed  that  the  difference  can  be  accounted  for  by  experimental 
error;  therefore,  the  results  are  encouraging.  Further  experimental  work  may 
be  done  using  a  search  coil  technique  to  measure  the  induced  current  across 
the  plate. 

FE  and  FFT  Methods  for  Pulsed  Magnetic  Fields.  As  discussed  in  the 
previous  section,  the  distribution  of  current  in  the  plate,  as  well  as  induced 
pressure  and  temperatures  are  for  harmonic  excitation.  For  pulsed  or  transient 
excitation  fields  B°(t),  one  may  decompose  the  field  into  its  spectral 


components 
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B  ( t )  e^dt 


(V-k) 


If  the  induced  current  is  calculated  as  a  function  of  frequency  J(o>),  for  an 
excitation  B0e-:iW^,  the  time  variation  of  current  can  be  found  from  the 


integral 


J(t)  =  f  J(M)  f(o))e"iWtda) 


(V-5) 


The  function  J(<jJ)  was  found  by  calculating  the  induced  current  for 
various  frequencies  using  the  nonlocal  theory  as  shown  in  Figure  V-ll. 

A  polynomial  was  then  fitted  to  these  points  over  the  significant  frequency 
domain  of  f(<Jd).  The  Fourier  integral  was  then  approximated  by  a  finite 
sum  and  the  summation  carried  out  using  a  fast  Fourier  transform  algorithm. 
The  resulting  time  history  of  current  in  the  plate  is  shown  in  Figure  V-10. 
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Figure  V-ll.  Spectrum  for  Current  at  x  *  -0.62A  for  Local  and 
Nonlocal  Solutions. 
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From  this  history  the  magnetic  force  and  induced  temperature  histories  in 
the  plate  can  be  calculated. 

Effect  of  Edges  on  Induced  Eddy  Currents.  One  of  the  immediate 
results  of  infrared  scan  experiments  on  induced  eddy  currents  in  plates  was 
the  observation  that  high  current  densities  will  occur  when  the  excitation 
field  is  concentrated  near  the  edges  of  the  plates.  This  can  be  seen  in 
Figure  V-l2 a  & b  where  a  coil  is  moved  closer  to  the  edge  of  a  rectangular 
plate.  The  results  of  the  one-dimensional  program  also  show  a  distinct  edge 
effect  for  both  the  paired  current  filament  induction.  Figure  V-9,  and  the 
uniform  field  induction.  Figure  V-13. 

Pressure  Distributions  Due  to  Tilted  Induction  Coils.  Both  magnetic 
forming  devices  and  magnetic  levitation  devices  use  current  carrying  coils 
near  plate-like  conductors.  The  effect  of  tilting  the  induction  coil  has 
been  calculated  as  shown  in  Figures  V-lU  and  V-15.  The  increased  magnetic 
pressure  under  the  filament  close  to  the  plate  produces  a  moment  on  the  plate 
and  a  restoring  moment  on  the  coil.  The  effect  of  lateral  movement  of  the 
coil  can  also  be  seen  in  Figure  V-15,  and  the  effect  of  lateral  movement  on 
the  restoring  moment  can  be  calculated. 

Effect  of  Reduced  Matrix  Band  on  Nonlocal  Solutions.  One  manifestation 

of  the  FE  implementation  of  the  nonlocal  theory  is  that  the  algebraic  FE 

equations  become  full  rather  than  banded.  Nevertheless,  the  nonlocal  effects 

2 

which  cause  this  loss  of  banding  are  proportional  to  1/r  in  which  r  is  the 
distance  from  a  self-field  source  point  on  the  sheet  conductor  to  the  field 
point  in  question,  eq.  (III-2).  One  method  of  restricting  the  nonlocal 
effect  at  any  field  point  on  the  sheet  would  be  to  exclude  source  points 
at  distances  from  the  field  point  greater  than  some  cut-off  value  of  r. 

For  a  uniformly  spaced  finite  element  mesh,  this  would  yield  a  banded  matrix. 


Figure  V-13.  Calculated  Edge  Effects  for  Uniform  Induction  Field. 
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Figure  V-lL.  Magnetic  Pressures  Due  to  a  Tilted  Pair  of  Current 
Filaments  Centered  Above  a  Long  Conducting  Plate 
(R  *  0.07,  tilt  *  10°). 


-56- 


-57- 


but  the  size  of  the  band  would  now  be  determined  by  the  cut-off  distance. 

The  choice  of  an  appropriate  cutoff  distance  will  be  dependent  upon  the 
particular  problem.  The  solution  to  the  restricted  nonlocal  problem  with 
an  appropriate  cutoff  distance  would  be  negligibly  different  from  the  full 
nonlocal  solution,  and  considerable  savings  in  computation  would  be  achieved 
due  to  a  reduction  in  the  nonlocal,  integrations  required,  eq.  (III-2). 

As  a  preliminary  test  of  this  strategy,  the  expedient  of  merely 
reducing  the  bandwidth  of  the  FE  equations  was  utilized  by  neglecting  all 
terms  outside  of  selected  bandwidths.  For  the  nonuniform  mesh  necessary  for 
the  eddy  currents  induced  by  a  pair  of  external  current  filaments,  this  pro¬ 
cedure  is  not  equivalent  to  the  selection  of  a  cutoff  distance.  (A  finer 
mesh  is  required  under  each  filament).  Nevertheless,  as  illustrated  in 
Figure  V-l6,  the  reduction  in  bandwidth  provides  a  variation  in  the  solution 
that  approaches  the  full  nonlocal  solution.  These  results  provide  evidence 
that  reduced  nonlocal  solutions  will  be  efficient  and  useful  in  some  situations. 
Future  research  will  follow  up  on  this  indication  by  developing  criteria  for 
the  selection  of  appropriate  cutoff  distances  and  by  implementing  this  strategy 
in  computer  programs. 

Two-Dimensional  Analysis 

As  described  in  Chapter  IV,  a  two-dimensional  finite  element  code  has 
been  written  to  predict  the  induced  eddy  currents  in  flat  plates  where  the 
skin  depth  is  equal  to  or  greater  than  the  plate  thickness.  This  code  has  been 
partially  verified  in  some  limiting  cases  such  as  the  low  magnetic  Reynolds 
number  limit.  In  this  limit  the  solution  for  the  stream  function  reduces 
to  a  Poisson  equation.  For  the  uniform  normal  external  magnetic  field,  the 
problem  is  analogous  to  the  torsion  of  a  shaft.  In  this  limit,  the  FEM  solution 
for  a  square  and  rectangular  plate  has  been  checked  with  analytical  solutions 
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for  torsion  of  rectangular  shafts. 

The  stream  function  contours  for  a  rectangular  plate  in  a  uniform 
magnetic  field  with  sinusoidal  time  variation  is  shown  in  Figure  V-17. 


The  resulting  comparison  of  induced  currents  across  the  middle  of  the  plate 
has  checked  the  one-dimensional  solution  in  the  low  R  limit  (Figure  V-l8). 

Nondestructive  Testing  (NDT)  Application.  One  of  the  immediate  appli¬ 
cations  of  the  two-dimensional  results  is  the  possibility  of  using  induced  eddy 
currents  and  infrared  scanning  to  detect  cracks  and  flaws  in  sheet-like  con¬ 
ductors  such  as  plates  and  shells .  When  the  current  flow  is  interrupted  by  a 
crack,  the  increased  current  concentration  near  the  crack  tips  increases  the 
induced  temperature  there  which  can  be  detected  by  an  infrared  scanner.  Such  an 
example  is  shown  in  Figure  V-19  where  the  tip  of  the  crack  or  slit  shows  up  as 
a  bright  "hot  spot". 

The  flow  of  current  around  a  crack  or  slit  in  a  plate  was  calculated 
using  the  two-dimensional  code  for  in  the  large  skin  depth  limit  (R  -*■  0)  for 
a  few  notch  widths.  The  stream  function  contours  for  one  case  are  shown  in 
Figure  V-20.  These  contours  show  the  flow  of  current  around  the  crack  tip. 
Contours  of  constant  temperature  are  also  shown  in  Figure  V-21.  These  results 
are  in  qualitative  agreement  with  the  infrared  measurements  since  they  clearly 
show  the  temperature  or  J“  hot  spots  near  the  edges  of  the  crack  or  slit. 

This  program  can  thus  be  used  to  develop  a  catalog  of  thermograms 
for  various  flaws,  cracks,  and  cutouts  in  sheet-like  conductors.  Also  the 
code  can  be  used  to  calibrate  eddy  current  probe  NDT  devices  by  predicting 
the  change  of  impedance  in  the  probe  as  the  eddy  current  probe  is  brought 


near  a  crack  or  flaw. 


Figure  V-17.  Stream  Function  Contours  for  a  Long  Rectangular 
Plate  Excited  by  a  Harmonic  Uniform  Field 
(R  =  0.0012). 
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Figure  V-lS.  Comparison  of  One-  and  Two-Dimensional  Solutions  of 
Currents  Induced  Across  the  Middle  of  a  Long 
Rectangular  Plate  by  a  Harmonic  Uniform  Field 

(R  =  0.03). 
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VI.  CONCLUSIONS 

This  technical  report  describes  the  accomplishments  of  the  first  year 
of  a  research  project  directed  toward  developing  numerical  methods  for  the 
coupled  analysis  of  forces,  currents,  and  stresses  in  thin  elastic  structures 
exposed  to  time-dependent  magnetic  fields.  A  stream  function  for  the  induced 
eddy  currents  is  used  to  reduce  the  three-dimensional  electromagnetic  problem 
to  an  integro-differential  equation  in  the  two-dimensional  surface  of  the  sheet 
conductor.  Although  formulations  are  presented  herein  for  both  flat  and 
cylindrically  curved  conductor  sheets,  to  date  the  research  has  focused  on 
the  response  of  flat  plates  to  steady  state  harmonic  and  pulsed  magnetic  fields. 

Finite  element  codes  for  one-  and  two-dimensional  flat  plates  have  been 
developed  which  can  calculate  the  induced  current  and  magnetic  fields  in  the 
plates,  the  magnetic  pressure  distribution  on  the  plates,  and  the  temperature 
distributions  in  the  plates  due  to  Joule  heating.  The  eddy  current  calculations 
include  the  nonlocal  or  self-field  effects.  These  programs  have  been  verified 
by  comparison  to  experimental  results  and  to  known  analytical  solutions. 

The  programs  have  been  applied  to  a  number  of  configurations.  The  re¬ 
sults  of  these  analyses  give  rise  to  the  following  conclusions: 

1.  The  numerical  method  based  on  the  stream  function  method  produces 
results  in  agreement  with  experimental  observations  and  analytical  solutions. 

The  success  of  this  method  is  encouraging  because  the  reduced  dimensionality 
promises  significantly  greater  efficiency  than  methods  based  on  a  full  three- 
dimensional  approach.  In  this  sense,  the  stream  function  method  is  analogous 
to  the  boundary  integral  equation  method  (boundary  element  method).  However, 
when  nonlocal  effects  are  taken  into  account  the  algebraic  finite  element  equa¬ 
tions  become  full  rather  than  banded,  and  this  increases  the  computational  cost. 
This  can  be  ameliorated  and  a  degree  of  bandedness  restored  by  a  proper  modi- 
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fication  of  the  nonlocal  integrals  to  account  for  the  rapid  decrease  of  the  self¬ 
field  effect  with  increasing  distance  between  sources  and  field  points. 

2.  For  medium  to  high  magnetic  Reynolds  numbers,  the  contributions  of  the 
self-field  effect  yield  significantly  different  eddy  currents  than  those  pre¬ 
dicted  by  a  purely  local  solution.  Most  previous  work  in  eddy  current  prediction 
and  magnetomechanics  has  neglected  this  effect.  Moreover,  for  a  harmonic  field, 
the  nonlocal  or  self-fields  effects  must  be  taken  into  account  to  obtain  the 
time-averaged  forces  on  the  conductor. 

3.  The  techniques  and  programs  that  have  been  developed  will  have  appli¬ 
cations  to  a  number  of  practical  problems.  One  such  application  illustrated  in 
this  report  is  the  use  of  induced  eddy  currents  for  the  nondestructive  testing 
(NDT)  for  cracks  and  flaws  in  sheet-like  conductors  such  as  plates  and  shells. 

The  programs  can  be  used  to  design  and  calibrate  eddy  current  probe  NDT  devices 
and  to  develop  a  catalog  of  thermograms  for  various  flaws.  Another  possible 
application  is  to  geomagnetic  prospecting,  in  which  very  large  induction  coils 
on  the  surface  of  the  earth  might  be  used  to  detect  the  presence,  location, 

and  orientation  of  seams  of  conducting  ores.  The  programs  could  be  used  to 
simulate  this  process  and  thereby  play  a  role  in  the  interpretation  of  prospecting 
data. 

The  present  computer  programs  can  be  used  to  predict  the  dynamic  forces 
on  flat  plates  when  the  motion  of  the  plate  is  small.  The  immediate  task  of 
future  research  is  to  extend  this  to  large  deformations  of  the  conductors  where 
true  coupling  of  the  magnetic  and  mechanical  behavior  occurs.  In  addition,  it 
is  planned  to  extend  this  work  to  curved  conducting  surfaces  such  as  cylindrical, 
toroidal,  and  other  shell  configurations.  Finally,  applications  to  nondestructive 
testing  are  being  explored  further. 
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APPENDIX  A 


Description  of  Eiqjerimental  Apparatus 


Description  of  UTI  Infrared  Scanner 

Commercially  available  infrared  detectors  include  point  sensing,  line, 
and  two-dimensional  or  plane  scanners.  Two  different  commercial  plane  scanners 
have  been  used  over  the  history  of  this  research  program,  though  for  the 
reported  phase  of  the  study  only  the  UTI  Scanner  was  used. 

The  scanner  consists  of  three  main  subsystems:  (l)  the  detector, 

(2)  the  continuous  grey  coding  of  temperatures  on  a  black  and  white  monitor, 
and  (3)  a  color  coding  and  color  display  monitor.  The  scanner  used  was  manu¬ 
factured  by  UTI  Corp.  of  Sunnydale,  California  based  on  an  earlier  design 
called  Spectrotherm  2000,  by  another  company.  The  detector  is  a  mercury- 
cadmium  telluride  (HgCdTe)  crystal  maintained  at  J8°K  by  liquid  nitrogen. 

HgCdTe  is  a  photoconduct ive  element  sensitive  to  infrared  radiation  in  the 
2-12  ym  wavelength  range.  Radiation  from  different  points  in  the  field  are 
focused  on  the  crystal  by  two  rotating  mirrors  for  both  horizontal  and  vertical 


The  sensitivity  of  this  unit  is  as  low  as  0.2°C.  The  output  from  the 
crystal  detector  is  used  to  produce  a  variable  grey  spot  on  a  cathode  ray 
display  tube.  White  to  black  grey  tones  can  be  used  for  hot  to  cold  or  visa 
versa.  An  isotherm  option  is  available  as  well  as  a  graphics  temperature 
display.  One  horizontal  scan  of  temperatures  is  graphed  below  the  two-dimensional 
grey  image  of  the  temperature  field. 

The  two-dimensional  grey  infrared  picture  has  525  raster  lines  and  600 
elements  per  line  on  a  7.6  cm  x  7.6  cm  screen.  The  field  can  be  focused  from 
10  cm  to  6  meters.  The  depth  of  field  is  15  cm  at  a  distance  of  80  cm. 
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The  scan  and  graphic  temperature  display  takes  2  seconds,  though  on  a 


never  model  a  one-second  option  is  available. 

The  electronic  data  can  be  stored  electronically  indefinitely  so  that 
pictures  can  be  taken  of  the  display  tube.  Also  a  zoom  capability  can  be 
used  to  magnify  and  examine  only  part  of  the  scanned  field. 

Finally  the  data  can  be  color  quantized  in  a  separate  unit.  The  UTI 
uses  a  Sony  Trinitron  (trade  name)  12.7  x  19-7  cm  display  screen  and  ten  colors. 
Each  color  represents  a  temperature  interval.  The  resulting  display  is  a 
beautiful  color  set  of  isotherms. 

The  AGA  Corp.  also  manufactures  a  two-dimensional  field  scanner 
system.  The  scan  time  is  shorter  (l6  per  second)  but  the  image  is  coarser 
(210  horizontal  lines,  lkO  elements  per  line)  and  no  storage  of  data  was  possible. 
This  unit  was  used  in  earlier  research,  however,  with  some  success. 

Specimen  Preparation 

The  ammount  of  radiation  transmitted  to  the  detector  depends  not  only 
on  the  temperature  of  the  transmitting  object  but  on  its  emissivity.  To  avoid 
spurious  data  due  to  differences  in  emissivity  the  conductors  were  sprayed  with 
a  high  emissivity,  black  coating.  Two  brands  used  were  NEXTEL  brand  $uede 
coating  (3101  series)  and  3M  Velvet  Coating-Optical  Black.  These  coatings 
also  decreased  spurious  infrared  reflections  from  other  objects. 

Search  Coil  Apparatus 

The  infrared  data  was  compared  with  experimentally  measured  eddy  currents 
by  using  a  differential  search  coil.  It  is  known  from  electromagnetics  that 
the  Jump  in  tangential  magnetic  field  ocross  a  thin  conductor  is  directly  pro¬ 
portional  to  the  current  in  the  conductor.  To  measure  the  tangential  magnetic 
field  at  the  surface  of  the  plate,  two  small  multi-turn  coils  were  wound  and 
aligned  with  their  axes  parallel  with  the  plate.  The  induced  currents  then 


generate  time-dependent  magnetic  fields  which  can  be  detected  by  these  coils . 
Since  the  output  voltage  is  proportional  to  the  change  in  magnetic  field,  a 
passive  R-C  integrating  circuit  was  used  in  series  with  the  coils. 

The  coils  were  connected  to  cancel  voltages.  Thus  the  primary  field 
of  the  induction  coil  would  produce  almost  zero  voltage  in  the  differential 
coils.  (Small  errors  due  to  the  finite  separation  of  the  coils  and  the  inhomo 
geneity  of  the  induction  coil  field  can  occur).  The  resulting  voltage  is 
almost  entirely  due  to  the  field  of  the  eddy  currents. 
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Technical  Paper 

K.Y.  Yuan,  F.C.  Moon,  and  J.F.  Abel,  "Magnetic  Forces  in  Plates  Using  Finite 
Elements",  Proceedings ,  Third  Engineering  Mechanics  Division  Specialty- 
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